Thermal conduction in molecular chains: Non-Markovian effects 
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We study the effect of non-Markovian reservoirs on the heat conduction properties of short to 
intermediate size molecular chains. Using classical molecular dynamics simulations, we show that the 
distance dependence of the heat current is determined not only by the molecular properties, rather 
it is also critically influenced by the spectral properties of the heat baths, for both harmonic and 
anharmonic molecular chains. For highly correlated reservoirs the current of an anharmonic chain 
may exceed the flux of the corresponding harmonic system. Our numerical results are accompanied 
by a simple single-mode heat conduction model that can capture the intricate distance dependence 
obtained numerically. 
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I. INTRODUCTION 

The problem of heat conduction through molecu- 
lar structures has recently attracted lot of attention 
[Hi S H, 0) [E [ll with potential applications in thermal 
machinery [7|, |8|, information processing and computa- 
tion [H, [lo| , and thermoelectricity [ll|, ■ 

One of the major open questions here is what are the 
factors that dominate thermal transport, the molecular 
structure, the contacts, or both [31? Another central issue 
is the determination of the system size [N) dependence 
of the heat current J. While Fourier's law suggests the 
relation ,/ cx N~^, extensive studies of heat flow in low 
dimensional systems have resolved a J cx N^*^ behav- 
ior, where a usually deviates from 1 p^ . Specifically, 
for harmonic chains one gets a = in the Markovian 
limit [l3|, i.e. the heat current does not depend on sys- 
tem size. This ballistic behavior results from the lack 
of scattering mechanisms between normal modes. In the 
harmonic limit the heat current thus reflects the spectral 
properties of the thermal reservoirs. In other words, it 
crucially depends on the details of the boundary condi- 
tions [i^, [la|. In contrast, in strongly anharmonic sys- 
tems where local thermal equilibrium exists, one expects 
that the steady state energy current will not depend on 
the properties of the contacts. 

Recent studies of heat transport in disordered low- 
dimensional harmonic chains have manifested the influ- 
ential role of the contacts' spectral properties on the 
asymptotic a value [l^ [l3] ■ Subsequent works exempli- 
fied this effect within anharmonic lattices [11 [H. While 
these works have typically focused on the asymptotic 
length behavior, our objective here is to systematically 
study the effect of the reservoirs' spectral properties on 
the thermal transport in short to intermediate size molec- 
ular junctions that are of experimental relevance. 

Using classical molecular dynamics simulations, we 
analyze the distance dependence of the current and 
the chain's temperature profile for Markovian and non- 
Markovian thermal baths considering either harmonic or 
anharmonic internal molecular interactions. We find that 



the spectral properties of the reservoirs play a crucial role 
in determining the size dependence of the thermal cur- 
rent. Thus, one should carefully interpret experimental 
results [3j] , as both molecular structure and the properties 
of the boundaries critically determine the junction con- 
ductivity. Another interesting finding is that for highly 
correlated noise, anharmonic chains conduct more effec- 
tively than the corresponding harmonic systems. We 
qualitatively explain our numerical results using a single- 
mode heat conduction model that can be solved analyti- 
cally [lS|2l|. 



II. MOLECULAR DYNAMICS SIMULATIONS 

We present here detailed classical molecular dynam- 
ics simulations of steady state heat transfer through 
one-dimensional (ID) molecular chains coupled to non- 
Markovian reservoirs. We model the molecule as a chain 
of N identical atoms. The end particles 1 and TV are 
connected to heat baths of temperatures Tl and Tr re- 
spectively. The dynamics is governed by the generalized 
Langevin equation 



Xk{t) = 



1 dHo 
m dxk 
1 OHq 
m dxi Jq 

° - / dt'jR{t^t')xNit') + mit)- 



, k = 2,3..N- 1 



- / dt'jLit-t')xi{t') + r]L{t), 



m dxN Jo 



(1) 



Xk is the position of the k particle of mass m, and pk [see 
Eq. (O] is the particle momentum. Hq is the internal 
molecular Hamiltonian. 7^, and 7^ are friction constants 
and rjL and 77^ are fluctuating forces that represent the 
effect of the thermal reservoirs. These terms are related 
through the fluctuation-dissipation relation {n = L, R) 
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where fc^ is the Boltzmann constant. We consider here an 
exponentially correlated Ornstein-Uhlenbeck (0-U) noise 



ln{t~t') 



\t-t'\/T; 



(3) 



with the intensity e and a correlation time Tc- For short 
correlation times the heat baths generate an uncorrelated 

(white) noise, 7„(t — t') — — ^ 2e„i5(t — t'). The Fourier 
transform of the 0-U correlation function, to be used 
below, is 



7n(w) 



2e„ 



1 + {LUT^y 



(4) 



A simple approach for implementing the 0-U noise in 
numerical simulations is to introduce auxiliary dynami- 
cal variabl es y ijt) and yN{t) for the L and R baths re- 
spectively [23 • The new equations of motion for the first 
particle are 
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where /iL(t) is a Gaussian white noise, (uLit)) = and 
{fiL{t)lJ-L{t')) = S{t — t'). An equivalent set of equations 
exists for the N particle, interacting with the R ther- 
mal bath. The coupled equations, Eq. ([1]) for particles 
2..N — 1 and ([5]) with its N equivalent, are integrated 
using the fourth order Runge-Kutta method to yield the 
positions and velocities of all particles. The heat flux can 
be calculated from the trajectory using [13] 



J = 



Af-l 



2(A^- 



{{vk+Vk+i)F{xk+i- Xk)), (6) 



We also use param eters where the anharmonic coefficient 
is large, a/^/mD of order 1. We refer to the later case 
as " anharmonic" . 

Unless otherwise stated, in the numerical simulation 
presented below we have typically used the following 
parameters: D — 367. S/i'^ kJ/mol, a = 1.875j^ 
Xeq = 1.54 A, and to— 12 gr/mole. These numbers de- 
scribe a c-c stretching mode for = 1 [iJ]. We take 
V — 0.01 for the harmonic model, while in the anhar- 
monic case we use v = Q. We also assume that the two 
reservoirs have the same type of spectral function (0-U) 
with equal strength e = e„ and noise correlation time 
Tc = tJ}. Depending on the situation, we have used inte- 
gration time step At = 10^^ — 10^* l/w, where oj is the 
molecular frequency in the harmonic limit. We also take 
care of the required inequality At <^ Tc- 
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FIG. 1; Distance dependence of the heat current in non- 
Mar kovian systems for harmonic (full), and anharmonic (dot- 
ted) models, (a) Gaussian white noise; (b) 0-U noise with 
Tc=8 xlO"^ ps; (c) 0-U noise with Tc = 0.01 ps; (d) 0-U 
noise with rc=0.Q4 ps. Tr = 3QQi«:, Tl ^OK,e = 50 ps"^ in 
all cases. 



where F{r) — —dHo{r)/dr, Vk = Pfc/w, and we average 
over time after steady state is achieved. 

We describe next the molecular structure of the chain. 
We model the interactions between the atoms using a 
Morse potential of dissociation energy D, width a, and 
an interatomic equilibrium separation Xeq, 



N 



k=l 



N 

2m ^ V 

k=l 



-a{xk + i-Xk-x^q) 



(7) 



We consider two sets of parameters: In the first case the 
potential width is taken to be very small a <C 1, so as 
the potential energy is practically harmonic with a force 
constant 2Da^. We refer to this model as "harmonic". 



Figure [T] presents the heat current for harmonic (Jh) 
and anharmonic (J a) chains calculated with different 
memory times Tc- Panel (a) shows the heat current in 
the Mar kovian limit. We find that the energy flux in har- 
monic systems does not depend on size, while it decays 
with distance for anharmonic chains in agreement with 
standard results [13]. When the noise correlation time is 
increased, an interesting behavior is observed [panel (b)] : 
While Jh remains a constant to a good approximation, 
the anharmonic flux manifests an initial rise, followed by 
a decay for long enough chains. As the memory time 
is further increased (c), the current of an anharmonic 
system saturates, and is approximately a constant over 
the relevant sizes. This observation interestingly shows a 
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FIG. 2: Resolving the distance dependence of the thermal 
current in non-Markovian O-U systems for harmonic (a) and 
anharmonic (b) models, Tc^OM ps, Tr = 300A', Tl = OK, 
e = 50 ps~\ 




FIG. 3: Decrease of heat current with increasing bath cor- 
relation time for harmonic systems. Tr — 300K, Tl = OK, 
e = 50 ps-^ N=2 (full), = 5 (dotted), = 10 (dashed). 
The inset zooms on the high Tc values. 



counteracting effect between the molecular contribution 
to the heat current and the reservoirs spectral proper- 
ties. For highly correlated reservoirs (d) both harmonic 
and anharmonic currents are slightly enhanced with dis- 
tance. Surprisingly, in this case the anharmonic junction 
conducts better than a fully harmonic system. 

We explain next these observations. First we clarify 
why Jh- and J a for short chains- increase with N for 
non-Markovian baths. As was shown in Ref. ^2§\ the 
dominant heat conducting vibrational modes of alkanc 
chains are shifted towards lower frequencies with increas- 
ing molecular size. Since within the O-U model "f{i^) 
(reflecting the system-bath coupling) is larger at lower 
frequencies, the current gets enhanced with distance. 

Next we explain the intricate current-distance behav- 
ior of anharmonic systems. Anharmonic interactions lead 
to scattering processes between the molecular modes. 
These scattering effects are more influential with increas- 
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FIG. 4: Decrease of heat current with increasing bath corre- 
lation time for an anharmonic molecular model. N=2 (full), 
iV = 5 (dotted), TV = 10 (dashed). Tr = 300 A", Tl = OK, 
e-n ~ 50 ps~^. The inset zooms on intermediate Tc values 
where J a is independent of length. 

ing chain length. In Markovian systems this results in 
the enhancement of the junction resistance, thus it leads 
to the reduction of current with N . However, in non- 
Markovian systems these scattering effects are actually 
beneficial for transferring energy from molecular modes 
which are above the reservoirs' cutoff frequencies, into 
low energy modes that overlap with the solids vibrations. 
The interplay between these two effects leads to a rich be- 
havior: If is higher than the molecular frequencies, 
here of the order of 150 ps^^, anharmonic effects lead 
to the reduction of current with size, see panels (a)-(b). 
In the opposite small cutoff limit (d), = 25 1/ps, 
harmonic systems can transfer only those modes that are 
in the reservoirs energy window, while anharmonic junc- 
tions better conduct by scattering high energy modes into 
low frequencies. For Tc ~ 0.01 ps the two effects practi- 
cally cancel and J a weakly depends on distance (c). 

Fig. [2] presents the distance dependence of the energy 
flux for both harmonic and anharmonic chains assuming 
solids of long memory time, Tc = 0.04 ps. While it is 
difficult to make a definite conclusion, we find that Jh 
and Ja obey different functional forms. 

Next we systematically explore the dependence of the 
heat flux on the reservoirs' correlation time. Figs. [3] 
and m manifest that the Markovian behavior sustains for 
times up to Tf, ~ 5 X 10""^ ps. For longer correlation 
times the heat current significantly decays with Tc for 
both harmonic and anharmonic chains. The inset of Fig. 
|4] further shows that for Tc ~ 6 x 10"'^ ps the anharmonic 
heat flux is practically distance independent, see also Fig. 
[Tfc). As discussed above, this intriguing behavior results 
from an effective cancellation between internal molecular 
interactions, leading to the decay of current with N, and 
the reservoirs properties, which can lead to an enhance- 
ment of current with distance. In the next section we 
present a simple analytical model that can capture this 
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intricate behavior. 

In Fig. [5] we study the temperature dependence of 
the heat current for both Markovian and non-Markovian 
chains for a representative length N=8. We find that 
both harmonic and anharmonic systems show a hnear 
current-temperature relationship in the range T = 0—300 
K. The thermal conductance {J/ AT) calculated from 
Figs. [l]|3]is therefore approximately independent of tem- 
perature. 

We have also analyzed in Fig. [6]the e dependence of the 
current for Markovian and 0-U systems, and found an 
approximate linear relation in the low dissipation regime. 
For Markovian baths the current decreases as when 
the coupling is strong, e > 100 ps~^. We expect that the 
colored noise model will demonstrate a similar behavior 
for very strong molecule-bath interactions [l9j . 

We conclude this section by noting that the qualita- 
tive behavior observed above (Figs. [l]|4]) applies for a 
broad range of temperatures and coupling parameters. 
We expect that similar characteristics will be discovered 
in molecular systems of various anharmonic internal in- 
teractions. 
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FIG. 6: Thermal flux as a function of system-bath coupling 
strength e, Tii=300 K, Tl=0 K. Parameters and lines setting 
are the same as in Fig. (5) 



quency luq. The total Hamiltonian includes three terms, 
H = Hb + Hq + Hj, where for a harmonic local mode 
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FIG. 5: Temperature dependence of the thermal flux, N=8, 
Tr—SOO K. Harmonic chain with white noise (full); Anhar- 
monic chain with white noise (dashed); Harmonic chain with 
O-U noise, rc=0.01 ps, e— 50 ps~^ (dashed-dotted); Anhar- 
monic chain with O-U noise, Tc— 0.01 ps, e— 50 ps~^ (dotted). 
The inset zooms on the O-U simulations. 



III. SINGLE-MODE HEAT CONDUCTION 
MODEL 

We present here a simple model that yields a qualita- 
tive explanation for the influential role of non-Markovian 
reservoirs on the heat transfer properties of short to 
intermediate {N — 2 — 20) molecular chains. In our 
simple picture, heat current in a linear molecular sys- 
tem is dominated by a specific vibrational mode of fre- 



Hp 



Hn 



2m., 

Pi 



H, = 



2M 2 



(8) 



Hb includes two thermal reservoirs n — L, R of different 
temperatures, each consisting a set of independent har- 
monic oscillators with masses mj, coordinates qj^n, and 
momenta pj^n ■ Hq represents the (single) relevant molec- 
ular mode with coordinate qo, momentum pq, mass M 
and frequency luq- It can be written equivalently in the 
energy representation as Hq = J2i=q i 2 ^"^olOi^li ^ = 1- 
If we sum over states up to infinity (Iloq Tn), this term 
describes a harmonic mode as in Eq. ([5]). We can also 
model an anharmonic molecule by truncating the sin- 
gle mode spectrum to include only few vibrational states 
pol . I21I . The system-bath interaction is taken to be bi- 
linear, with Xj^n as the coupling constant. 

Assuming weak molecule-bath coupling at both ends, 
an analytical expression for the heat current in steady 
state can be derived using the master equation formalism 
[20I [2H . In the harmonic limit and at high temperatures 
(T > luq) the thermal current is given by a Landauer type 
expression [2^ 



J 7L(^0)7fl(^0) , 

•JH = — 7 — ^- — tKb AJ , 

where AT = Tr-Tl and 7„(cj) = f E 



(9) 



is the Fourier transform of the friction constant 7„(i) 
[26| . The memory damping can be also expressed 
in terms of the reservoir's spectral function gn{^) = 



5 



-S{uj — ujj) as 7„(w) — g„(tj)/4ijj. Equation 



Mrrijujj ^ 

(121) clearly manifests that in the harmonic limit the heat 
current is exclusively determined by the spectral proper- 
ties of the reservoirs. This behavior results from the lack 
of scattering mechanisms in the harmonic system. 

As mentioned above, within this simple picture we can 
also model a highly anharmonic molecule by truncating 
the single mode spectrum. For a two-level model (/=0,1) 
at high temperatures the heat current is given by [20| 



7i('^o)7ij(wo) ^0 



AT, 



(10) 



7l(wo) + 7-r(wo) 2Tb 

where the temperature of the local model (sometimes re- 
ferred to as a bridge B) is 

7L(wo)rL + 1r{i~oq)Tr 



Tb = 



(11) 



7fl(wo) 
TU)) provide the heat cur- 



Though expressions © and 
rent for a simplistic model, they may serve us for gaining 
qualitative understanding of heat transfer in an A^-sites 
molecule. The key element here is the observation that 
in short linear chains relatively few modes contribute to 
the transport of thermal energy [2^. For simplicity, we 
may assume that a single mode dominates the dynamics, 
and use the following generic form to describe its size 
dependence 



LOq w ujm 1 + 



N 



(12) 



lom is the asymptotic frequency for large TV, /? is a con- 
stant that is specific for the material. Note that Eq. (fT^ 
does not necessarily describe the variation of the low- 
est vibrational mode of the chain with increasing length. 
Rather, this is a qualitative expression for the variation of 
the dominant heat conducting mode with size: For short 
chains the spectrum is significantly altered with size. For 
large enough chains wq is approximately fixed. 

We can clearly discern in Eqs. (|5|l- (ll0p the role of dif- 
ferent factors (contacts and internal interactions) in de- 
termining the heat current. While the reservoirs spectral 
properties enter these expressions through the damping 
rate 7, calculated at the relevant molecular frequency 
ujQ, anharmonic interactions yield an extra ujq/Tb fac- 
tor that accounts for the local mode occupancy. Since 
both of these terms depend on frequency, thus on size 
through Eq. (jUl), the resulting N dependence is not 
trivial. We study next various models for the damp- 
ing element, Markovian and non-Mar kovian, assuming 
for simplicity that the system is symmetric with respect 
to the two ends, 7 — 7„ {n = L, R). 

I. Markovian process. For a Markovian process the 
spectral function is Ohmic, g{uj) = Aeuj, thus the friction 
is a constant, 7(0;) = e. In this Hmit Eqs. ^ and pU)) 
reduce to 
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FIG. 7: Qualitative behavior of the heat current in the 
single-mode heat conduction model for Markovian and non- 
Markovian reservoirs, (a)-(b) Heat current for a Gaussian 
white noise [Eq. (I13|l ]: (c)-(d) Current for an O-U noise with 
a short memory time, Eq. (I14|l with ujmTc < 1; (e)-(f) Cur- 
rent for an O-U noise with a long memory time, Eq. (I14p with 

iOMTc ~ 1. 
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FIG. 8: Qualitative behavior of heat current in the single- 
mode heat conduction model with a power law spectral func- 
tion IlSl), s = 1. 



These expressions are consistent with standard results 
[l3| : The current in harmonic systems does not depend 
on size, while in anharmonic models it decays with dis- 
tance. In Fig. [7] (a)-(b) we qualitatively exemplify this 
behavior by simulating Eq. (jl3p , reproducing the results 
of panels (a)-(b) of Fig. [1] 

II. Ornstein-Uhlenbeck noise. The exponentially cor- 
related (O-U) process "f{t) — e/Tce~'*'/'^'= transforms into 
a Lorentzian in frequency domain, 7(0-1) = The 
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fluxes ([9]) and (fTO|) then become (e = e„ and Tc = t") 



IV. TEMPERATURE PROFILE 
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where we derived the approximate distance dependence 
by assuming that woTc > 1 and taking /3/iV < 1. These 
expressions clearly demonstrate an enhancement of the 
current with N for both harmonic and anharmonic sys- 
tems when the reservoirs have (each) a highly correlated 
noise. 

In Fig. [7] (c)-(f) we simulate Eq. HI]). For short 
correlation time we reproduce the results of Fig. [Ijb), 
manifesting an enhancement of the anharmonic flux fol- 
lowed by a decrease of current. For very long correlation 
times the current systematically increases with size, in 
agreement with the numerical data, Fig. [1] (d) . 

///. Power law models. The spectral function g{uj) — 



j/u. 



is a widely accepted description of solids, 



euj e 

is the reservoir cutoff frequency. Unlike the 0-U noise 
which satisfies the differential equations (O, this noise 
cannot be reduced into a multi-component Markovian 
process, thus it is not trivial to simulate \2f\. Within the 
single-mode heat conduction model the heat current for 
harmonic and anharmonic modes ([5|l- pl}|) is given by 



J 



J 



H 

(P) 



AT 



'2{Tl + Tr) 



(15) 



where loq is the dominant frequency for heat transport 
(fT2)) . For an Ohmic {s ~ \) bath, taking u)c ^ ojq, 
the dynamics reduces to the Markovian limit (fT5)) , where 

J\j is weakly enhanced with distance and J\ ' decays 
with iV. This behavior is exemplified in Fig. [8] (a)-(b). 
In contrast, for ~ "-^Oj a different trend is observed: 
while monotonically increases with length (c), in 

anharmonic systems the current first increases with size, 
then falls down (d). For s = 2 we expect that, quite 
interestingly, the heat current in harmonic models will 
decay with distance. This is because for long chains the 
dominant conducting frequencies are red shifted, while 
the solid spectral function is maximal at Wc- 

We conclude that short to intermediate size molecu- 
lar chains coupled to general environments can mani- 
fest a rich behavior: The thermal current can either in- 
crease or decrease with size, critically depending both on 
the molecular internal interactions and on the reservoirs 
spectral properties. Despite its simplicity, the model pre- 
sented here provides a useful starting point for explaining 
the complicated behavior observed within atomistic clas- 
sical molecular dynamics simulations. In the next section 
we show that this model is also a useful tool for estimat- 
ing the temperature profile along molecule. 



We investigate next the effect of the reservoirs spectral 
properties on the local temperature profile in harmonic 
and anharmonic chains. Within the single mode heat 
conduction model one can define in steady state the tem- 
perature of the local mode as /csTs ~ c^o X); ^^ii where 
Pi is the population of the / vibrational state. In the clas- 
sical limit it can be shown that this expression reduces 
to Eq. (HB, Tb = 7i.(^Q)T^+7H(^o)Tjt ^^j. ^^^^j^ harmonic 

and anharmonic local modes [20,, i21|] . We therefore ex- 
pect that when the reservoirs have the same type of spec- 
tral function with the same coefficients, the temperature 
profile is independent of the type, and is given by the 
arithmetic average Tq = {Ti^+Tfj) /2. On the other hand, 
when the two reservoirs have different types of spectral 
densities (or different memory times), the temperature 
profile along the chain depends on the model. For exam- 
ple, if Tl < Tr and 7h < 7l, Tb ^ IbTrJ-^i^ < Tr. 
In contrast, for ^ Tr but > 7^, Tb ^ Tr. This 
behavior was observed by Saito et al. in the harmonic 
limit using a quantum master equation formalism (T6| . 

Using classical molecular dynamics simulations we 
compute the temperature profile by calculating the mean 

kinetic energy for each particle, = (^)jr-- Here is 
the momentum of the fcth particle calculatecf after steady 
state develops, and we average over time and initial con- 
figurations using the procedure and parameters of section 
II. 

Fig. [5] manifests that the bulk temperature profile is 
close to the arithmetic average (150 K) for both Marko- 
vian and non-Markovian systems when the reservoirs are 
identical. At the boundaries we see interesting features 
noted earlier in Refs. [ij, [l^, [H, [2^: The tempera- 
ture close to the cold end iV) slightly rises, even to val- 
ues higher than the average temperature. For an an- 
harmonic system a small temperature gradient develops 
at the chain center. This effect is more pronounced for 
Markovian reservoirs. Note that in Markovian harmonic 
systems the temperature at the chain center, e.g. Tio, is 
higher than the anharmonic value. This trend is reversed 
for 0-U baths, as here anharmonic interactions facilitate 
equilibration with the slow environment. 

We employ next reservoirs of different spectral proper- 
ties, and study the heat current and temperature profile 
in this asymmetric system. We assume that both reser- 
voirs are of 0-U type, but take different memory times 
7^ ■''(f^- Fig. [in] manifests that the current in a har- 
monic system is independent of the direction of the asym- 
metry, i.e. if the cold bath is weakly or strongly coupled 
to the molecule. In contrast, the current in an anhar- 
monic chain is significantly altered when the asymmetry 
is reversed: It is higher when the hot bath is slow. This 
response to spatial asymmetry is the underlying principle 
of the "thermal rectifier" 0, [1^. 

Fig. [TT] shows that the bridge temperature dramati- 
cally responds to the direction of the asymmetry for hoth 



FIG. 9: Temperature profile along harmonic (full) and an- 
harmonic (dashed) chains coupled to reservoirs of identical 
spectral properties, (a) Gaussian white noise; (b) O-U noise 
with Tc = 0.04 ps. N=20, Tr^SOO K, Tl=0 K and e = 50 
ps~^ in all cases. 
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FIG. 10: Heat current in molecular chains coupled to reser- 
voirs of different spectral properties 7^ t^, el = en =50 
ps"^ ri^ = 0.002 ps, ri* = 0.04 ps: harmonic (full); an- 
harmonic (dashed), — 0.04 ps, = 0.002 ps: harmonic 
(dashed-dotted); anharmonic (dotted). rjf,=300 K and Tl=G 
K in all cases. 

harmonic and anharmonic systems. We find that it is 
close to the temperature of the reservoir with the short 
memory time. In other words, the molecule better equi- 
librates with the Markovian reservoir. This behavior is 
consistent with the results of Ref . [l^ , and qualitatively 
agrees with Eq. pdj). 



V. DISCUSSION AND SUMMARY 

In this paper we have investigated the effect of the 
contacts' spectral properties on the thermal conduction 
of harmonic and anharmonic short to intermediate size 
chains using classical molecular dynamics simulations. 

When the system size is smaller than the mean free 
path, heat conduction is dominated by harmonic inter- 
actions. In such systems the spectral properties of the 



280^ ^ ^ ^ ^ 

5 10 15 20 

k 

FIG. 11: Temperature profile along harmonic (full) and an- 
harmonic (dashed) chains coupled to reservoirs of different 
spectral properties 7^ r^, el = =50 ps^^. = 0.04 
ps, ri- = 0.002 ps (top); = 0.002 ps, = 0.04 ps (bot- 
tom). The temperature profile was calculated for an A'^ = 20 
sites chain using rfl=300 K, Tl=G K in all cases. 



contacts play a central role in determining the dynamics 
in the junction. We found that the distance dependence 
of the current strongly varies with the baths (noise) cor- 
relation time: While for Markovian baths the current in 
harmonic systems is independent of size, it slightly in- 
creases with length when the reservoirs have long mem- 
ory. The effect is even more dramatic in anharmonic 
chains. Here the current decays with N in the Marko- 
vian limit, while it gets enhanced with distance when the 
reservoirs have long memory times. 

Other interesting observations are: (i) We identified a 
regime where the anharmonic flux is practically distance 
independent due to the counteracting effects of molecular 
nonlinearities and bath correlations, (ii) For long mem- 
ory times the current through anharmonic chains might 
be higher than the analogous harmonic current. Thus, 
anharmonic interactions might play a surprising role, en- 
hancing the heat flux across molecules coupled to non- 
Markovian reservoirs. We also showed that the single- 
mode heat conduction model developed in Ref. [20| can 
essentially capture the length dependence and the tem- 
perature profile obtained within classical molecular dy- 
namics simulations. 

Before we conclude we briefly discuss the relationship 
of our calculations to relevant experiments. Wang et 
al. have recently measured the thermal conductance 
(/C = J/^AT, A is the cross section) of 8-10 alkandithiols 
monolayers Q. The values K. ^ 25-28 MWm^^K^^ were 
obtained, roughly independent of system size. Previous 
calculations 25] confirm this observation: The heat cur- 
rent of alkane chains is expected to be distance indepen- 
dent, even for very long (iV=100) molecules, due to weak 
anharmonic internal interactions in the system. More 
quantitatively, Eq. ^ predicts thermal conductance of 
/C 280 MWm~^K~^, using Markovian thermal baths 
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with e=10 ps""'^ and a cross section ^ = (5 x 10~^'')^ni~^. 
As discussed above, memory effects can reduce this num- 
ber by 1-2 orders of magnitude. Therefore the conduc- 
tance calculated here is in a reasonable agreement with 
experimental results. 

Several future directions are of interest. First, we 
have restricted our discussion to the 0-U correlation 
function, since it can be easily simulated as a multi- 
component Markovian process [23| . It is worth extending 
this study to include more realistic environments of uj'^ 
spectral properties. This type of noise correlations can 
be simulated using the inverse Fourier transform tech- 
nique [13, Hll, [13] ■ Other interesting questions are what 
is the role of the contacts in higher dimensions, and how 



do quantum mechanical effects alter the phenomena ad- 
dressed in this paper [H, M, M, M, M,M,E^- 

Understanding thermal transport in molecule-solid in- 
terfaces is important for molecular electronic applications 
[3^ . It is also an imperative step in the endeavor for de- 
veloping unique molecular level thermal devices [sl. [ol, IToj . 
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